Study of phase stability of MnCr using the 
augmented space recursion based orbital 
peeling technique 



Rudra Banerjee b Abhijit Mookerjee 



b,a 



^Department of Materials Science, S N Bose National Centre for Basic 
Sciences, Block JD, Sector III, Salt Lake, Kolkata 700098, India 

b Advanced Materials Research Unit, S N Bose National Centre for Basic 
Sciences, Block JD, Sector III, Salt Lake, Kolkata 700098, India 



Abstract 

In an earlier communication we have developed a recursion based approach to the 
study of phase stability and transition of binary alloys [1]. We had combined the re- 
cursion method introduced by Haydock, Heine and Kelly [2] and the our augmented 
space approach [3] with the orbital peeling technique proposed by Burke [4] to de- 
termine the small energy differences involved in the discussion of phase stability. 
We extend that methodology for the study of MnCr alloys. 
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1 Introduction 



In an earlier communication [1] we had introduced the augmented space re- 
cursion (ASR) [3] based orbital peeling method (OP) [4] as an useful and 
numerically accurate method for the calculation of the 'pair energies'. This 
allowed us to map the binary alloy problem onto an effective Ising-like model 
and study the stability or otherwise of different ordered phases that might 
arise if the disordered alloy is cooled below some critical temperature. The 
aim of this paper is to extend the use of this method to study MnCr alloys. 

These pair energies are small differences of relatively large electronic energies 
and a brute force calculation is likely to yield these small numbers with errors 
which are usually of the order of or larger than the numbers themselves. Direct 
estimations like the OP method are therefore appropriate in these situations. 
Our earlier analysis of FeCr, FePd and PdV alloys illustrated the success 
of the ASR-OP method in predicting the low temperature phases of these 
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Fig. 1. (Color Online) Densities of states for a series of alloy compositions for 50-50 
MnCr alloy. 

binary alloys and their stability. The analysis involved ordering and mixing 
energies obtained from a generalized perturbation framework (GPM) [5,6] and 

the Fourier transform V(k) of the pair energies which are related to diffuse 
scattering intensities as described by Krivoglatz, Clapp and Moss (KCM) [7]. 
The basic ideas behind these methods have been described in detail in the 
references quoted above. For the sake of completeness, we shall briefly dwell 
upon the main points used in our analysis. 

The basis of our subsequent analysis is the electronic structure calculation on 
the MnCr alloy system. We have chosen the tight-binding linear muffin-tin 
orbitals (TB-LMTO) approach [8] which provides a first-principles, density 
functional based tight-binding sparse representation of the Hamiltonian. Such 
a Hamiltonian is appropriate for describing the electronic structure of random 
substitutional alloys in which the configurational fluctuations due to disorder 
are local. The sparseness of the TB-LMTO Hamiltonian is a suitable input for 
the recursion method. 

Realistic models for disordered alloys require us to go beyond the single-site 
mean-field approximations and treat effects of configuration fluctuations more 
accurately. This we shall attempt through the ASR [9,10]. As mentioned ear- 
lier, the pair-energies have to be accurately determined. We shall adopt the 
recursion based OP for such calculations. The approach followed here will 
provide a unified recursion based methodology to address such problems. 

In this communication our choice of system is the MnCr alloy. We shall first 
obtain the configuration averaged density of states for a series of alloy compo- 
sitions using the TB-LMTO-ASR. The ASR has been discussed in great detail 
earlier. We shall only stress here that it generalizes the single-site mean-field 
approximations to include the effect of the configuration fluctuations of the 
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near neighbourhood of a site and yields configuration averaged Green func- 
tions which retain their analytic and lattice translational symmetries even 
after approximation. The results for the 50-50 alloy are shown in Fig. 1. The 
converged potential parameters are input into the OP calculations that follow. 

We shall organize the paper as follows : in section 2 we shall describe the OP- 
ASR method for obtaining the effective pair energies, the expression for the 
ordering and mixing energies and the analysis for the Fourier transform of the 
pair energies. In the section 3 we shall analyze the results for the equi-atomic 
MnCr alloy. Concluding remarks will be given in the final section. 



2 Methodology 



2. 1 Total energy and pair energies 



The simplest model which analyzes the emergence of long ranged order from 
a disordered phase is the Ising model. Our approach will attempt to map 
the energetics of the binary alloy problem onto an equivalent 'spin-half Ising 
model. We need a derivation of the lowest configurational energy for the alloy 
system in terms of effective multi-site interactions, in particular "effective pair 
energies" (EPE) [1 1] . We need to accurately and reliably determine the EPE. 
Our approach will be to start with the disordered phase, set up a perturbation 
in the form of concentration fluctuations associated with an ordered phase and, 
from its response, study whether the alloy can sustain such a perturbation. 
This approach includes the generalized perturbation method (GPM) [12], the 
embedded cluster method (ECM) [13] and the concentration wave approach 
[14]. 

We shall begin with a homogeneously disordered alloy A x Bi_ x , where every 
site is occupied by either an A or a B type of atom with probabilities pro- 
portional to their concentrations. We define the 'occupation' variable to 
be a random variable which takes on the values 1 and according to whether 

the site labeled R is occupied by an A or a B atom. Its average <Cn^;^>= x. 
This perturbative approach expands the total internal energy of a particular 
atomic configuration as follows : 

E = V^+Y^V^ K? + ^E E V M 6n R Sn Rf + ... (1) 

R R R'ytR 

here 8ng = — x and <^8n^>— 0. If the configuration is homogeneously 

disordered then it immediately follows that <CE^>= E<n s = V^. From the 
above definition we can interpret the other two expansion terms as follows : 
if E 1 is the configuration averaged total energy of a configuration in which 
any arbitrary site labeled R is occupied by a atom of the type I and the 
other sites are randomly occupied, and E IJ is the averaged total energy of 

— * — * 

another configuration in which the sites R and R' are occupied by atoms of 
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the types I and J respectively and all other sites are randomly occupied, then 
from equation (1) it follows that : 



V ( } ] = E A — E B Vi 2 l = E BB + E AA - E AB - E BA (2) 

R RR 1 

The one-site energy is unimportant for bulk ordered structures emerging 
from disorder. It is important for emergence of inhomogeneous disorder at 



r(2) 

a v 

factors governing emergence of bulk ordering. 



surfaces and interfaces [15]. The pair energies V~l, are the most important 



The interpretation of equation (2) immediately allows us to introduce a method 
to obtain the pair potentials directly rather than calculate the total ener- 
gies and then subtract them. Since they are small differences (of the order of 
mRyd) of large energies (of the order of 10 3 Ryd), a direct calculation will pro- 
duce errors larger than the differences themselves. The orbital peeling method 
(OP) [4] based on recursion [2] was introduced by Burke precisely to calculate 
such small differences, albeit in a different situation. 

The total energy of a solid may be separated into two terms : a one-electron 
band contribution Ebs and an electrostatic term Ees which includes several 
contributions : the Coulomb repulsion of the ion cores, the correction for dou- 
ble counting terms due to electron-electron interaction in E BS and a Madelung 
energy in case the model of the alloy has atomic spheres which are not charge 
neutral. The renormalized cluster interactions defined in equation (1) should, 
in principle, include both Ebs and Ees contributions. Since the renormalized 
cluster interactions involve the difference of cluster energies, it is usually as- 
sumed that the electrostatic terms cancel out and only the band contribution 
is important. Obviously, such an assumption is not rigorously true, but it has 
been shown to be approximately valid in a number of alloy systems [6]. We 
shall accept such an assumption and our stability arguments starting from the 
disordered side, will be based on the band structure contribution alone. 

The effective pair interactions can be related to the change in the configuration 
averaged local density of states : 

V M = J dE (E - E F ) An(E) (3) 



where An(E) is given by : 

i AB 

An(E) = 3mV Tr <£/ - H^)' 1 ^ 

71 u 



£u = 25 u — 1, i.e. is ± 1 according to whether I — J or I ^ J. There are 
four possible pairs IJ : AA, AB, BA and BB. H^ J ^ is the Hamiltonian of 
a system where all sites except R and R' are randomly occupied. The sites 

— * — * 

labeled R and R' are occupied by atoms of the type / and J. This change in 
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the averaged local density of states can be related to the generalized phase 
shift 77(E) through the equation : 



An(E) 



dr](E) 
dE 



d_ f det < G AA {E) > det < G BB (E) > 
dE \ ° g det < G AB (E) » det < G BA (E) > 



G 7J (£') is the resolvent of the Hamiltonian H( /J ). The generalized phase shift 
i](E) can be calculated following the orbital peeling method of Burke [4]. We 
shall quote only the final result : The pair energy function is defined as : 



AB L 



f M ,(E) = f(R - R>) = £ E iu j dE ' ( E> - E ) log « ^ J (£') » 



IJ a=l 



^max 

EE 

7J a=l 



E Z fe a ' /J - E Pk IJ + (Np' U - N a / J ) E 
k=i 



k=i 



(4) 



<C G 7J ^> denote the configuration averaged resolvents in which the 
orbitals from L = 1 to (a — 1) are deleted. Z%' and P%' are its zeros 
and poles and N^' IJ and ATp' IJ are the number of such zeros and poles of 
< G 7J (^) > below The zeros and poles are obtained directly from the 
recursion coefficients for the averaged resolvents and these are obtained from 
the TB-LMTO-ASR. This method of zeros and poles enables one to carry 
out the integration in equation (3) easily avoiding the multi-valuedness of the 
integrand involved in the evaluation of the integral by parts. 

The pair energy is then given by : V { ?1 = fa#(E F ) = V (2) {R - R!). The last 
expression follows if the background system, in which the A and B type atoms 

— * — * 

are immersed at R and R', is homogeneously disordered. 



2.2 The Krivoglatz-Clapp-Moss analysis 



Philhours and Hall [17] have suggested, and Clapp and Moss [7] have formally 
shown that a sufficient (but not necessary) condition for a stable ground state 
is that the wave vector of concentration waves corresponding to an ordered 
phase lie in the positions of the minima of the Fourier transform of the pair 
energy function 

V{k)= E exp {ik ■ {R - R')} V^(R-R') 

R-R' 



The above statement follows from the expression for the inverse susceptibility 
which measures the response of the disordered system to the concentration 
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fluctuation perturbation described above. 
X-^ocl + ^l-x^KfrW 



In a zero-th approximation V e g(k) = V(k). Corrections to the effective pair- 
function has been described in detail by Chepulskii and Bugaev [18]. We have 
used here the Ring Approximation suggested by the authors as the one most 
suitable for our analysis : 

V cS {k) = V(k) - 09/2)(l - 2xf J F(q)F(k-q) (5) 



where 

V(q) 



F(q) 



l + x(l-x)pV(q) 
2.3 Ordering and mixing energies 

Finally, the GPM expression also gives the ordering energy : 

AE ord = - VonQn 



where n is a n-th nearest neighbour of an arbitrarily chosen site (which we 
label 0) and Q n = (x/2)(N BB - xN n ), N BB is the number of BB pairs and 
N n the total number of pairs in the n-th nearest neighbour shell of 0. 

With reference to the total energies of the pure constituents, in the approxi- 
mation where we only restrict ourselves to pair energies and ignore all three 
body energies and higher, the so-called mixing energy is given by : 

&E mtx = -\x{\ - x)J2 N n v£ ] 



The averaging procedure ASR has been described in great detail in many 
earlier papers and we refer the reader to the review [9] in which the method and 
its relation to the CPA and its generalizations has been discussed extensively. 



3 Results and Discussion 



We have calculated the composition dependent pair energy functions using 
the TB-LMTO-ASR coupled with the orbital peeling technique and the result 
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Fig. 2. (Color Online) Pair functions calculated for MnsoCrso, variation shown 
against E-E^, calculated from TB-LMTO-ASR-OP. (top) Nearest neighbour pair 
function at a distance y/3a/2, a being the equilibrium lattice parameter (bottom) 
second, third and fourth nearest neighbour pair functions at distances a, \[2a and 
\/3a 
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Fig. 3. Pair energies calculated for MnCr at a 50-50 composition. 

for the equi-atomic composition is shown in Fig. 2 The nearest neighbour pair 

energy function fi(R — R', E) shows the characteristic shape of a positive lobe, 
indicating ordering, near the position of half filling, flanked by negative lobes 
indicating segregation near empty and complete filling fractions. 



We should note that in our approach, both the pair energy function itself and 
the position of the Fermi energy depend upon the composition of the alloy 
and its band filling. This is in contrast to some analysis (like the Connolly- 
Williams) which depend on similar, but composition independent, pair energy 

functions. The Fig. 3 shows the pair energies V^ 1 {\R — R'\) for the equi-atomic 
composition. 
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Table 1 

Weights for different neighbouring shells for seven different bcc based equi-atomic 
superstructures. 



Neighbouring shells 



Structure 


1 
1 
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q 
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Segregated 


1 nn 

x.uu 


n IK 
U. ( O 


l.uU 


q nn 


B2 


-1.00 


0.75 


1.50 


-3.00 


B32 


0.00 


-0.75 


1.50 


0.00 


Bll 


0.00 


0.25 


-0.50 


0.00 


ST1 


-0.50 


0.25 


0..00 


0.50 


ST2 


0.00 


-0.25 


-0.50 


0.00 


ST3 


0.50 


0.25 


-0.00 


-0.50 



The behaviour of the pair energies for MnCr indicate ordering tendencies 
up to the fourth nearest neighbours. The pair energies rapidly converge to 
zero with distance. In fact, although we had calculated the pair energies up 
to the seventh nearest neighbour shells, their values beyond the fourth shell 
were smaller than the error bars of our calculational method and therefore 
these numbers were not really reliable and were not used for our analysis. 
The same is true for the results of the ordering energies for seven different 
structures and super-structures based on the body-centered cubic lattice. We 
have, therefore, calculated the ordering energies with contributions only up 



0.15 



E 0.1 
o 



>. 0.05 

a. 
E 

> 0.0 



lu -0.05 
m 



■a -0.1 
O 



-0.15 



SG 



B11 



B32 



j &2 



ST3 

1 



ST1 



B2 



Fig. 4. (Color Online) Ordering energies for seven different structures for MnCr 
based on Table 1. 
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Fig. 5. (Color Online) V e s(k) and contours on the plane bounded by (top) (001) 
and (100) and (bottom) (001) and (110) 

to the fourth neighbouring shell. The Table 1 gives us the weights Q n for the 
seven bcc based structures required to obtain the ordering energies in this alloy 
system. These superstructures are described in detail by Finel and Ducastelle 
[19]. 

Based on the Table 1 we present the ordering energies for the same seven 
bcc based structures for MnCr in Fig. 4. Unlike our earlier study of FeCr 
which indicated segregation and possible ordering in the ST3 superstructure, 
for MnCr ordering in the B2 structure with a possible competition from the 
superstructure ST1 is the most energetically probable. The contrast between 
the two alloys FeCr and MnCr is interesting : Fe segregates with Cr while Mn 
orders. The stainless steel alloys, a class of which are ternary alloys FeCrMn 
should then see competition between segregation and ordering. We intend to 
study the ternary compositions subsequently. 

Finally we shall examine the Fourier transform V(k) of the pair energies and 
carry out a Clapp-Moss type of analysis. The contour diagrams for Fourier 

transform V(k) in the is plotted in Fig. 5 for the (001) and (110) planes. 

The minimum occurs at k — (001) This is indicative of a possible B2 type of 
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ordering, as indicated by the ordering energy analysis. 

In order to ascertain whether the ordered state associated with the minimum 
is stable compared to the segregated species, we have to carry out a much 
more detailed analysis including the contribution of the energy of mixing. The 
mixing energy for the B2 structure is 0.16412 mRyd/atom This indicates that 
ordering is stable against segregation. 

We should comment here that the above discussion is based only on the elec- 
tronic contribution to the diffuse scattering intensity. At the temperatures that 
we are interested in there are contributions from the vibrational excitations 
in the system. Our aim here was to indicate the possibility of ordering ten- 
dencies in MnCr rather than accurate estimation of energetics and transition 
temperatures. In any detailed and accurate statistical mechanical calculations 
we must include the contribution of vibrational excitations to the free energy 
of the alloy. 



4 Conclusion 



In a previous work [1] we had introduced and examined the suitability and ac- 
curacy of the Augmented Space Recursion (ASR) based Orbital Peeling (OP) 
method for the generation of pair energies. In this communication we have 
extended these ideas into body-centered cubic MnCr alloy. We have looked 
at the phase stability of the equi-atomic alloy. Unlike our previous work on 
FeCr, FePd and PtV alloys, for MnCr we had no earlier theoretical work to 
compare with. However, since our applications to FeCr, FePd and PtV gave us 
satisfactory results, in good agreement with experimental evidence, we have 
confidence in our present results. With FeCr, this work will form the back- 
ground of our extension of this work to ternary FeMnCr stainless steel alloys. 
The recursion and ASR is now available with both relativistic corrections (in- 
cluding spin-orbit terms) [20] and non-collinear magnetism [21,22]. Since earlier 
we have proposed the ASR as a analyticity preserving generalization of the 
single-site mean-field coherent potential approaches, this work will provide 
further incentive to extend the use of the ASR to problems beyond the simple 
density of states calculations and in problems were relativistic corrections and 
non-collinear magnetism will play significant roles. 
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